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ABSTRACT 

For Very Long Baseline Interferometry (VLBI), the fringe spacing is extremely narrow compared to the field of view imposed by the 
primary beam of each element. This means that an extremely large number of resolution units can potentially be imaged from a single 
observation. 

We implement and test a technique for efficiently and accurately imaging large VLBI datasets. The DiFX software correlator is used 
to generate a dataset with extremely high time and frequency resolution. This large dataset is then transformed and averaged multiple 
times to generate many smaller datasets, each with a phase centre located at a different area of interest. 

Results of an 8.4 GHz four-station VLBI observation of a field containing multiple sources are presented. Observations of the cali- 
brator 3C345 were used for preliminary tests of accuracy of the shifting algorithm. A high level of accuracy was achieved, making 
the method suitable even for the most demanding astrometric VLBI observations. One target source (1320+299 A) was detected and 
was used as a phase-reference calibrator in searching for further detections. An image containing 13 billion pixels was constructed by 
independently imaging 782 visibility datasets covering the entire primary beam of the array. 
Current implementations of this algorithm and possible future developments in VLBI data analysis are discussed. 



Key words. Instrumentation: interferometers - Methods: data analysis - 
- Techniques: interferometric - Techniques: image processing 
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1. Introduction 

For an interferometer, the fringe spacing or angular resolution 
6 is approximately A/D where A is the wavelength and D is the 
baseline length. The field of view is dictated by the solid angle of 
the sky to which the antenna is sensitive. For parabolic antennas 
this can be approximated by a symmetric Gaussian with a half- 
power beam width of the primary beam given by 1 .2A/d where 
d is the diameter of the parabola. For VLBI, D/d and therefore 
the number of resolution units across the primary beam 0/(9 is 
particularly high (See tableQ}. For imaging purposes the number 
of pixels will be ~3 times higher in both dimensions than the 
number of resolution units to ensure proper sampling of the dirty 
beam. 

This means that very large images can be generated. 
However there are three caveats: (1) this technique is only use- 
ful if there are several sources sufficiently bright and compact for 
VLBI detection within the primary beam; (2) it is necessary to 
correlate the data with sufficient resolution in time and frequency 
to avoid bandwidth and time-average smearing; (3) there must 
be sufficient computing resources and appropriate techniques to 
handle the resulting large datasets. 

Progress in all of these three areas has been made in 
the decade since post-processing computers have been able 



Table 1 : Comparison of field of view in resolution units for dif- 
ferent arrays. 



Array 


d(m) 


D(km) 


D/d 


VLA" 


25 


36 


1440 


Merlin* 


32 


217 


8 680 


EfMaMcNo r 


100 


1334 


10334 


EVN 


100 


10180 


101 800 


VLBA 


25 


8611 


344440 



Notes. Diameter of largest antenna d; length of longest baseline D; 
the approximate number of resolution units across the primary beam 
D/d « 0/0. 

(a) A configuration. <A) Without the Lovell 76-m. {c> Array described in 
Sect.l3~n 
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to handle the maxi mum output rates of hardware correlators. 
iGarrett et al.l (11999) review early efforts at wide-field imaging 
of objects including gravitational len ses and multiple supernova 
remnants located in the same galaxy dPedlar et al.ll 19991) . 

Hardware correlators generally place some constraint on the 
number of output channels and on the output integration tim e 
(For example the JIVE Mark 4 correlator: see Camp beilll2004b . 
This in turn limits the possible field of view of the output dataset 
due to time and bandwidth smearing (see Sect. 12. 41 . Multiple 
sources can be imaged by using multiple passes through the cor- 
relator dRioia et al. 1 19991) however this is infeasible for a large 
number of sources. 
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In spite of these limitations further progress was made 
including VLBI o bservations of the Hubble deep field 
dGarrett et all 1200 ll) and the NOAO bootes field dGarrett et all 
2005). Le nc et alT (l2008) took data from low-frequency observa- 
tions and analysed it using wide-field VLBI techniques (a mix- 
ture of multiple correlator passes and UV shifting). This led to an 
unprecedented 28 detections. Nonetheless, time and frequency 
resolution constraints necessitated the discarding of data from 
longer baselines with increasing distance from the phase cen- 
tre; thus none of the above studies were able to image the entire 
area permitted by the field of view of the interferometer elements 
with full sensitivity and UV coverage. 

Constraints on the time or frequency resolution are much les s 
stringe nt in software corr elators such as the DiFX software cor- 
relator dDeller et al.ll2007l ) and are thus a good match to the re- 
quirements of wide-field imaging. The first attempt to image ob- 
jects acro ss the entire primary bea m without discarding data was 
made bvfMiddelber g et al.l (120081) . 

Technical developments in VLBI mean that bandwidths con- 
tinue to increase. The LBA currently works at a data rate 
of 512 Mbit s" 1 with higher rates possible for some antennas 
(Phill ips etai]l2009h. The EVN is p ushing towards 4Gbits~ 1 
from 1 Gbit s" 1 (Ivan Langveldd2009l) . The ongoing VLB A sen- 



sitivity upgrade will increase the bandwidth fr om 5 12 Mbit s to 
2 Gbit s -1 and then to 4 Gbit s" 1 dWalker et al.ll2007l) . increasing 
continuum sensitivity by a factor of approximately three. This 
will increase the density of detectable sources on the sky. 



2. Methods 

2.1. Correlation and generation of a visibility dataset 

In VLBI correlation it is necessary to generate a delay model 
to align the incoming sampled data as if the antennas were all 
on a plane perpend icular to the point on the sky being observed 
dSovers et al.l ll998). In order to achieve the accuracy required 
to ensure picosecond timing precision (required for microarc- 
second precision astrometry) the model used is necessarily ex- 
tremely complex. 

The delay model dictates the 'phase centre' of the data: the 
location where the residual fringe delay and rate are zero. If the 
field only contains a single point source at the phase centre, the 
visibility phases will be static along the time and frequency axes 
(ignoring un-modelled instrumental and atmospheric effects). A 
single source that is far from this phase centre will have visibili- 
ties that rotate rapidly in phase along the frequency and time axes 
due to the residual delay and residual delay rate respectively. 

If the data are over-averaged this will lead to a reduction 
in flux density, and increasing uncertainty in the position of the 
source. In an FX correlator, fixed blocks of samples are selected 
for each baseline and Fourier-transformed before correlation. 
Only samples in corresponding FFT blocks can be compared 
leading to fur ther reduced se nsitivity to sources distant from the 
phase centre (iRomnevll 1999b . 

The baseline vectors (u, v, w) are calculated along with de- 
lays. These define the baseline with respect to the source for 
each time integration and baseline and are used when imaging to 
compute the residual visibility phase A(f> for sources away from 
the phase centre as 



2tt 

A(f> — —(lu + mv) 



Traditionally u and v have been calculated using a purely ge- 
ometric model. This neglects many effects, the most significant 
of which is differential aberration, a relativistic effect due to the 
orbital velocity of the Earth (v©) with a magnitude of ~ v ffi /c, up 
to about 10~ 4 , implying image registration errors of ~ 1 pixel for 
shifts of 10 4 pixels if not corrected. This will lead to significant 
errors for wide-field imaging. 

More accurate values of u and v can be derived using the 
delay model. These are calculated by sampling the delay at two 
further points: one with a small offset in Z, another with a small 
offset in m. Along with the delay already calculated for the phase 
centre, the change in delay with respect to / and m can be calcu- 
lated to give baseline vectors which make equation (Q} accurate 
to first order and consistent with the delay model used: 



Idr dr 

( "' v ' w) = c| ¥'^' T 



(2) 



(1) 



for a source displaced (I, m) radians from the phase centre as a 
substitute for computing the full delay at that point. 



An additional benefit of using equation (0 to calculate base- 
line vectors is that it provides a well-defined prescription even 
in the case of very near-field targets, such as planets, asteroids 
and spacecraft. The DiFX software correlator produces data with 
improved baseline vectors as described above. 

The output of the correlator is a dataset consisting of (1) the 
right ascension and declination of the phase centre; (2) baseline 
vectors for each baseline and time integration; (3) the visibil- 
ity data themselves. These visibility data are 4-dimensional with 
baseline, time, frequency and polarisation axes. Each datum is a 
complex number representing phase and amplitude. Other data 
usually included are the delay model used and various data use- 
ful for calibration and imaging. 

2.2. Wide-field VLBI data reduction 

An obvious approach for imaging o ver the entire beam is to use 
the large visibility dataset directl y (lGarrettetaDl200ll) at least 
for an initial search for detections dLenc et al.l2008l) . Data reduc- 
tion algorithms such as TMAGR' and 'VTESS' in AIPS are able 
to account properly for time and bandwidth smearing, working 
on datasets with multiple channels, and even shifting the phase 
centre of the image . Some efforts are being made to parallelise 
these algorithms dBhatnagarl l2009) however no currently avail- 
able software suitable for VLBI data reduction allows this as 
far as we are aware. Even with the fastest CPUs currently avail- 
able data processing time becomes infeasible for VLBI images 
larger than a small fraction of the primary beam. Moreover, di- 
rect wide-field imaging imposes added complications such as 
correcting smearing due to non-coplanar effects which quickly 
lead to overwhelming computing requirements. Finally as shown 
in Sect. 12.11 current datasets provide information on how the de- 
lay varies across the wide-field image only to first order. This 
will inevitably lead to errors in position for sources far from the 
phase centre unless a sophisticated delay model is used for cor- 
rection at the imaging stage. 

Provided that the baseband data has been recorded (not nec- 
essarily the case for e-VLBI) it is possible to correlate more than 
once using a different delay model each time. This permits the 
generation of several datasets phase-centred on different areas of 
the primary beam. However, unless the correlator is capable of 
correlating multiple phase centres simultaneously, the correlator 
time required will soon become prohibitive for a large number 
of widely-separated sources. 
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2.3. Transforming visibility datasets 

Another possibility is to correlate once with time and frequency 
resolution sufficient to avoid significant smearing on any part of 
the field that is to be imaged. This large single dataset can then be 
transformed to be phase-centred on a different point. Since time- 
average and bandwidth smearing effects increase radially from 
the phase centre, the data can be averaged in time and frequency 
before writing to disk. Similarly, non-coplanar effects also in- 
crease radially from the phase centre and so are mitigated. This 
shifting can be done multiple times, in parallel if necessary, gen- 
erating multiple datasets. These can then be imaged separately. 



2.3.1 . The transform 

In order to transform the data to a new phase centre we must: 
(1) replace the phase centre coordinates; (2) replace the baseline 
vectors with those calculated for the new phase centre; (3) add 
a phase shift to every visibility. The shift is time, frequency and 
baseline dependent but the same for all polarisations. 

The dataset can then be averaged in time and frequency. The 
degree of averaging possible will depend on the image size de- 
sired. If all of the primary beam is to be imaged, the averaging 
will depend on the spacing of the phase centres. If individual 
sources are being picked out, it will depend on the source struc- 
ture and the uncertainty in the a priori position of the source. 
In any case uncertainties in the antenna delay calibrations typ- 
ically limit averaging to ~5 seconds in time and ~0.5 MHz in 
frequency for most centimetre wavelength observations. 

2.3.2. Limitations of current software 

Interferometry data reduction packages commonly have algo- 
rithms for performing this shift, for example AIPS has the task 
'UVFIX'. However these algorithms are designed for shorter 
shifts and positional accuracy is lost on longer shifts. This is 
dem onstrated in Sect . 14.21 and has also been noted by others 
(see iLenc et all 120081 Midderberg et al.1 120081) . The averaging 
must be performed separately from shifting and the data must 
be sorted on disk multiple times requiring more disk usage. 
Moreover, extremely large datasets cannot be handled by AIPS 
efficiently, if at alfl 

2.3.3. Accurate UV shifting 

By generating two delay models, one for the correlation phase 
centre and another for the desired phase centre we can transform 
the data very accurately. 

Consider n antennas located at different time-varying posi- 
tions with respect a source at phase centre s. Each antenna re- 
ceives and records a time-varying voltage v n (t). During corre- 
lation the delay models T„(f) will be used to shift the functions 
v„(t) to v„(? + t„(0) so that a wavefront emanating from s will 
be aligned for all n. We want to take a dataset correlated with a 
phase centre at s' (with delay models r'^f)) and transform it so 
that it is identical to a dataset correlated with phase centre s. 

Consider a baseline of antennas i and j. Points in v,-(f) and 
Vj(t) which would have been adjacent after the application of t„ 
are separated by St after the application of T'j(t) and T^(f)- This 
separation St is the therefore the delay that we will have to apply 
in order to transform our data to the new phase centre. To a first 



approximation it is given by Stq which is simply the difference 
between the two delay models for any time t : 



«T = (^(0-T»(0)-(T , /0-T/0). 



(3) 



However this simplistic correction assumes that the delay 
models are not changing with time (i.e. t„ — f' n — where 
the dot notation is used for time derivatives). Essentially, we 
have estimated the delay correction at time t, when in fact we 
should have estimated the correction for each antenna n at time 
t + (j'Jf) - r„(f)) in order to take into account the change in the 
delay model in the interval St. Thus, approximating the delay 
model as a linear function over this short span of time, a better 
approximation of St can be made: 

*r, = (7#) - r,(f))d + t;(f)) - (r'/t) - r 7 (r))(l + , (4) 

where we have also assumed that f„ = r'„. For ground-based 
VLBI, 6to and 6t\ typically differ by approximately 1 part in 
10 6 , while the difference between St\ and St is less than 1 part 
in 10 8 . The latter is shown empirically in Sect. 14.21 

It is quite common to refine the delay model in post- 
processing using data not available at the correlator at correla- 
tion time (examples include using improved Earth orientation 
parameters). All delay calibration should properly account for 
the time-variation of the delay model. However in most cases 
the shifts in the delay model are small enough that equation (0 
is sufficient. 

Finally the delay must be applied to the (Fourier trans- 
formed) visibilities. The phase shift <p to be applied is simply 

4> = 2n(y ■ St mod 1) (5) 

for each sky frequency v. 

2.4. Correcting for amplitude losses 

Amplitude loss to due to time-average and bandwidth smear- 
ing were de scribed as averag i ng los se s in Sect. 12. 1 1 an d else- 
where (see iBridle & Schwab! (1 19991) : iThompson etail {2001, 
Chap. 6.4)). 

The time-average smearing of each shifted visibility ampli- 
tude (R s ) compared to the amplitude of the same source at the 
phase centre (Ro) is given by 



Ro 



sine 



2tt- 



T ■ ti, 



(6) 



where f,„, is the pre-shift accumulation period and t is the time 
derivative of St. 

In an FX correlator such as DiFX where only a single fast 
Fourier transform ( FFT) is performed on e ach sample, the trian- 
gular lag function (Thompson et al. 2001, Fig. 8.17) results in a 
reduced visibility amplitude for an offset continuum source. This 
is given by 



— = 1 - St ■ Sv 

Ro 



(7) 



where Sv is the width of each spectral channefl The lag function 
entirely dictates bandwidth smearing effects for an FX correla- 
tor. 



1 Specifically there is a limitation on n IF x n channe i s 
rently 32768) when loading the dataset. 



'■polarisations 



(cur- 



2 Assuming the baseband data is sampled at the Nyquist rate and an 
FFT size in samples of twice the number of spectral channels 
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We can also approximate the bandwidth smearing by anal- 
ogy with the time-average smearing case (i.e. assuming a rect- 
angular bandpass for each spectral channel). Assuming that 
Ro - R s <K 1 and replacing f • f,„, in equation (O with St ■ 5v 
we get 



l 



(7T • 6t ■ 6vf 



(8) 



This is far less than that predicted by equation (|7} and makes a 
case for overlapping FFTs in an FX correlator. 

For full accountability, both delay models should be stored 
for each visibility in the output dataset along with the pre- 
shifted integration time and spectral channel width and the FFT 
setup of the correlator. This would allow full correction in post- 
processing. In the results we present in Sects. 14.31 and 14.41 all 
amplitude losses have been kept sufficiently small to be ignored. 

2.5. Correcting for the primary beam response 

Imaging to the edge of the primary beam is done routinely with 
connected-element arrays. In this case a correction is applied to 
the intensity across the final image. However for heterogeneous 
arrays it is more appropriate to apply an antenna by antenna cor- 
rection. This is achievable using the multiple phase centre tech- 
nique if we assume that the primary beam response is constant 
across each sub-image. The primary beam power function of a 
baseline is given by 



A 12 (l, m) = F*M, m) x F 2 (l, m) 



(9) 



where F\ and F2 are the voltage patterns of the two antennas in 
the m plan e (this is discussed in the context of VLBI arrays in 
IStrorrEOQl . 

The primary beam corrected amplitude of each visibility I,p 
on the baseline i, j in a dataset phase-centred in the I, m plane 
can be calculated from the observed amplitude I t, s : 



hid, m) 



*-obs 



AM, m) 



(10) 



where the amplitude of the power function has been scaled to 1 
at the centre of the primary beam. 

It is also possible to calculate the theoretical RMS noise 
limit: 



A7 m = 



1 



7s 



V2 • Av • T h 



\Aij(l,m)\ 2 



fi=N-l,j=N 

^ , SEFD, x SEFD: 



-1/2 



(11) 



In Section 14.41 a simple truncated sine function is used to 
model the primary beam. The limitations of this simple model 
should be noted. First, the primary beam is different at each 
frequency. Second, the primary beam correction may be dif- 
ferent for each polarisation due to offset feeds. This, and any 
other non-radial function for the voltage, will lead to a time- 
varying primary beam for antennas with an alt-azimuth mount. 
A full treatment of thes e effects for the VLBA is discussed in 
lMiddelbergetaD(l20Tol) . 

3. Observations of the 1320+299 Complex 

With the aim of testing the performance of the UV shifting 
method described above, we organised a four-station ad hoc 
VLBI observation. The target field selected was that of the 
quasar 1320+299 (4C 29.48). The data recorded at each station 



were processed with the DiFX software correlator installed at 
the Istituto di Radioastronomia (IRA), Bologna, Italy. Below we 
describe the target source, the observation parameters and the 
initial correlation. 

This source was first observed at 1 .4 G Hz and 5 GHz using 
the Westerbork Synthesis Radio Telescope dFeretti et al.l 119821) 
and later with the V LA in A and B confi gurations at 1 .4 GHz, 
5 GHz and 15 GHz dCornwell et al.l[l98 6. herein referred to as 
Cornwel l et al.). The sour ce was also observed by VERA at 
22 GHz dPetrov et al.ll2007l) with the component B chosen as the 
phase centre, however no detection was reported. 

VLA and WSRT observations reveal three aligned compo- 
nents: a 20th magnitude quasar is coincident with component A 
(see Fig. 1, in Comwell et al.). FigureQ] which we have gener- 
ated from data in the VLA archive shows the relative locations 
of the three components. 

Component A has a one-sided morphology: a prominent ra- 
dio core, unresolved by the highest-resolution observations car- 
ried out so far (Cornwell et al. VLA A configuration at 15 GHz). 
It has peak brightness 210mJy beam 1 at 5 GHz and an inte- 
grated spectral index a of -0.42 (where S oc V). A weak outer 
component is connected to the core region by a jet-like exten- 
sion. 

Component B has an asymmetric structure with a 'jet' which 
bends towards the north-east, extending to a distance of approx- 
imately 3" from the peak. It has a peak brightness of 75 mJy 
beam 1 at 5 GHz and is also unresolved by the VLA in A con- 
figuration. The integrated spectral index of this component is 
-1.12. 

The VLA A-configuration 5 GHz image of component C re- 
veals it to have a feature which extends in the north-west direc- 
tion. It has a peak brightness of 24mJy beam 1 at 5 GHz and an 
integrated spectral index of -0.92. 

The three components are polarised and the degree of polar- 
isation at maximum brightness ranges between 5 and 19%. No 
significant depolarization or rotation measure is seen. 

As a whole, the structure of 1320+299 looks rather peculiar. 
There is no conclusive evidence for the three components A, B 
and C to be physically related; none of the images from prior ob- 
servations show any bridge of emission connecting the compo- 
nents. On the other hand, according to Cornwell et al. the prob- 
ability of an unrelated source with A2Q cm flux density between 
250 and 450 mJy lying within ~25" of the central component B 
is only ~5x 10 . Cornwell et al. do not reach any firm conclu- 
sion about the classification of the components of the 1320+299 
complex. Our observation, the first dedicated VLBI observation 
of 1320+299 to our knowledge, is also aimed at a better under- 
standing of this puzzling complex. 

3.1. Observation 

We observed the 1320+299 complex making use of an ad hoc 
network of antennas, namely Effelsberg (100-m), Medicina (32- 
m), Matera (20-m), and Wettzell (20-m). The data were recorded 
using four base-band converters of the MarklV acquisition sys- 
tem and MK5A disk recorder with a sample rate of 16 x 10 6 
samples s -1 and four-level sampling for all baseband channels. 
The four observing frequencies were 8405 .49, 84 1 3 .49, 842 1 .49, 
and 8429.49 MHz. Only the upper side band of 8 MHz was 
recorded for each channel. The observing session started on 12th 
December 2007 at 05:00:00 UT and stopped at 08:00:00 UT. We 
were limited in our observing time by having only 1 Terabyte of 
storage available at the IRA computing centre at the time. 
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03.0 02.5 02.0 01.5 

RIGHT ASCENSION (J2000) 



Fig. 1: Image of the 1320+299 showing the three components. 
VLA in A configuration at 4.8 GHz; synthesised beam shown at 
the bottom left; contour levels -0.7, 0.7, 1.4, 3, 8, 20, 40, 100, 
200 mJy beam -1 . 



The observations were carried out using the phase-reference 
technique: switching between the target source 1320+299 and 
the calibrator with a duty cycle of 3 minutes on target and one 
minute on a point source phase calibrator. The phase reference 
technique was required to increase our chances to detect fringes 
from the individual components A, B, and C and as a possible 
first epoch for future measurements of absolute component mo- 
tions. The VLBA calibrator J1329+3154 (flux density 0.91 Jy at 
8.4 GHz), approximately 2.5 °from the 1320+299 complex, was 
used as phase-reference source. 

3.2. Correlation with the DiFX software correlator at IRA 

The MK5A disk packs recorded at each station were sent to the 
Medicina Radio Observatory, read back on a MK5A unit and 
loaded onto hard disks. These were then transferred to IRA. 

The correlation of the data was carried out on a 12-node clus- 
ter. Each node consists of two quad-core CPUs and the nodes 
are connected by 10 Gigabit Ethernet. Connected to the cluster 
is a RAID system served to the cluster using the GPFS file sys- 
tem. Throughout the correlation the CPU and network utilisation 
was very low since the bottleneck was reading the data from the 
RAID. With the modest observing parameters used for this ex- 
periment we were able to correlate in approximately real time. 



4. Analysis 

4.1. Implementation of the shifting algorithm 

CALC 9.1 is typically used to generate delay models for DiFX. 
The algorithm described in Sect. I2.1l was implemented by mod- 
ifying the DiFX FITS conversion tool to allow a second delay 
model to be read. This could be used to apply any kind of modi- 
fied delay model to a correlated dataset with high precision. 

4.2. Verification of the shifting algorithm 

First it was necessary to test the accuracy of the UV shifting 
technique described in sect. 12.3.31 The most meaningful measure 
is to compare a dataset correlated on phase centre s and shifted 



to phase centre s' with a dataset correlated directly with phase 
centre s' . 

To compare the tw o resulting datasets the method outlined 
in Tingav et al.l (120091) for comparing two correlators was used: 
a single sub-band is extracted and the differences between the 
phases of the two datasets are measured and assessed for any 
systematic error that could degrade the astrometric precision 
even after calibration and averaging. This analysis is carried out 
twice: first the two datasets are vector-averaged along the time 
axis before comparison, then the process is repeated with spec- 
tral averaging. In each case the shifted dataset is subtracted from 
the reference dataset. The error is the mean value of this new 
dataset. 

A scan of the source 3C345 was used for this test. Imaging 
this source showed it to be unresolved with a flux density of 5 
Jy. One minute of data were correlated twice: once on the known 
position of the source and once with a 1' offset added to both the 
R.A. and Dec coordinates (equivalent to ~ 10 000 turns of phase 
for each of the three baselines). The second dataset was then 
shifted back, using both the AIPS task UVFIX and our algo- 
rithm. Further tests were then done with shifts 1 and 2 orders of 
magnitude larger respectively. The smallest of the three shifts is 
similar in magnitude to those required to obtain the results pre- 
sented in Sect. 14.41 The largest is equivalent to shifting beyond 
the edge of the primary beam for 10 m dishes on a 10000 km 
baseline. In each case the correlation parameters were chosen to 
keep time-average and bandwidth smearing to a level of approx- 
imately 5% for the first two shifts and 30% for the final shiffl 

All of the datasets were averaged post-shifting to a 2 s inte- 
gration time and 32 channels per sub-band. The outer 8 channels 
and first and last time integrations were discarded before analy- 
sis. 

Figures [2] and [3] illustrate the results of these tests for a sin- 
gle baseline for the shorter shift. The results for all shifts are 
summarised in table [2] The r 1 value is given by 



(12) 



where 6<p is the residual phase error between the two datasets 
and / is the expected value of 6(p from the fit. 

The low r 2 value for the linear regression suggests that the 
linear regression does not explain the variation of the points; in 
other words, the residual error is mostly random. Moreover there 
is no sign of any systematic error with increasing shift. The small 
random error is consistent with the expected noise introduced by 
smearing and would be an insignificant error even in the most 
precise VLBI measurements. 

Table [3] compares the difference in amplitudes compared 
with those predicted in equations [6] and Q Since the source is 
very strong, the large amplitude loss in the final shift still leaves 
a large SNR and so there is still good agreement between the 
phases in the two datasets. 

4.3. Searching for the three components 

A full correlation of the target data was made with the phase cen- 
tre close to the a priori coordinates of B (the pointing centre). 
An integration time of 0.2 s and 512 channels per sub-band kept 
time-average and bandwidth smearing to an extremely low level. 



3 The actual correlation parameters were: 128, 1024 and 2048 chan- 
nels per 8 MHz sub-band and integration times of 0.4, 0.032 and 
0.0064 s for the three shifts. 
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Fig. 2: Comparison of a dataset shifted using a second correlator delay model (our algorithm) with a UVFIX-shifted dataset (a) and 
with the reference dataset (b) for a single baseline (Matera-Medicina) with a positional shift of 1' N and 1' W. All spectral channels 
have been vector averaged. In (a) crosses denote the UVFIX data and circles the dataset shifted with our algorithm. In (b) crosses 
denote the dataset shifted with our algorithm and circles the reference dataset. 
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Fig. 3: Comparison of a dataset shifted using a second correlator delay model (our algorithm) with a UVFIX-shifted dataset (a) and 
with the reference dataset (b) for a single baseline (Matera-Medicina) with a positional shift of 1' N and 1' W. All time integrations 
have been vector averaged. In (a) crosses denote the UVFIX data and circles dataset shifted with our algorithm. In (b) crosses denote 
the dataset shifted with our algorithm and circles the reference dataset. 



Table 2: Accuracy of phase shifts. 



Baseline 




Shift 1 








Shift 2 








Shift 3 








Shift/turns 


Error/ 


r 2 


2 
r. 

lime 


Shift/turns 


Error/ 


r 2 


r 2 

tune 


Shift/turns 


Error/ 


r 2 

freq 


2 
r. 

time 


MaMc 


5.5 x W 


-0.18 


0.017 


0.087 


5.5 x 10 4 


-0.21 


0.034 


0.077 


5.4 x 10 5 


1.90 


0.46 


0.408 


MaWz 


11 x 10 3 


0.29 


0.044 


0.008 


1 1 x 10 4 


-0.76 


0.086 


0.150 


11 x 10 5 


-0.59 


0.13 


0.007 


McWz 


5.5 x 10 3 


-0.08 


0.021 


0.017 


5.6 x 10 4 


0.04 


0.002 


0.051 


5.6 x 10 5 


-0.97 


0.25 


0.242 



Notes. Shift magnitudes (v-6t), errors and r 2 values, for 3 different shifts done using a second delay model. The r 2 is calculated from a least-squares 
fit to a linear model after vector averaging each spectral channel in time (r 2 , ) and after averaging all spectral channels in each time integration 



The correction given in equation (Ql was deemed unnecessary. 
Three shifted datasets were then generated centred on the a pri- 
ori coordinates of the three components. Imaging of the phase- 
reference calibrator showed it to be entirely unresolved. The 
three shifted datasets were then calibrated in the standard way 



using the phase reference source and then imaged. Only compo- 
nent A was detected and after primary beam correction and self- 
calibration was revealed to be entirely unresolved (see Fig. [U 
with a flux density of 103 mJy beam -1 (VLA observations in 
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Table 3: Accuracy of amplitude correction. 



Baseline 




Shift 1 






Shift 2 






Shift 3 






time 


lag 


total 


actual 


time 


lag 


total 


actual 


time 


lag 


total 


actual 


MaMc 


1.63 


4.1 


5.7 


5.4 


1.05 


5.1 


6.1 


5.7 


4.44 


25.2 


28.5 


28.8 


MaWz 


0.42 


8.2 


8.6 


7.9 


0.28 


10.3 


10.6 


10.5 


1.38 


51.4 


52.1 


55.0 


McWz 


0.39 


4.1 


4.5 


4.4 


0.25 


5.2 


5.4 


5.2 


0.89 


26.2 


26.9 


28.6 



Notes. Percentage amplitude errors for 3 different shifts predicted due to (1) time-average smearing; (2) the triangular lag function of an FX 
correlator; (3) the product of (1) and (2). (4) The actual loss in amplitude observed averaged across all spectral points and time integrations. 




13 23 00.877 00.876 



Fig. 4: Map of 1320+299A centred on RA 13 h 23 m s :8733 
Decl. +29°4r44'.'814 after primary beam correction and self- 
calibration, imaged with a large circular restoring beam 
(7.5x7.5 mas); peak flux 102 mJy beam -1 ; contour levels -4, 4, 
8, 16, 24, 32, 48, 64, 80, 96mJy beam- 1 . 



Comwell et al. imply a flux density of approximately 170mJy 
beairT 1 ). No detection was made of the other two sources. 



4.4. Imaging the entire primary beam 

Finally, in order to search for any sources not present in previous 
studies over the relatively wide areas over which flux was recov- 
ered by Comwell et al.and to test the feasibility of using UV 
shifting to generate extremely large images, the entire primary 
beam was imaged. The 'edge' of the primary beam for the array 
was arbitrarily chosen to be the point at which the Ef-Mc base- 
line had the same predicted sensitivity as the Wz-Ma baseline 
- their relative sensitivities as a function of offset being calcu- 
lated using the summation in equation (fTTT i. assuming a radial 
sine function for each antenna's voltage pattern. This leads to a 
radius of 90" from the pointing centre, approximately double the 
area enclosed by the HPB W of Effelsberg. 

Shifting in turn to phase centres placed on a grid, and imag- 
ing each one allows the imaging of this entire field with each 
sub-image being produced from a small dataset. Spacing the 
phase centres 5" apart in declination and 0.5 s (equivalent to 
6.5") apart in right ascension gives a convenient grid and results 
in a total of 782 datasets (each 20MB in size) required to image 



the entire primary beam. With images of 4096 x4096 pixels with 
each pixel having dimensions of 0.85 x 0.7 mas, the entire area 
is covered with a slight overlap between the images. Averaging 
each sub-band to a single channel incurs bandwidth smearing of 
4% on this image. Losses due to non-coplanar effects are ~1%. 

In view of the strong detection of component A, two re- 
finements were made to the calibration procedure: (1) phase- 
reference solutions were derived from a fringe fit on component 
aQ; (2) an attempt was made to subtract a simple model of com- 
ponent A from the UV datasets; this was found to reduce spuri- 
ous bright sources significantly, particularly in the fields closest 
to component A. 

The calibration a nd imaging of each sub-field was automated 
using ParselTongue dKettenis et alJ2006h . Since much of the cal- 
ibration could be reused for each field, this process proceeded 
very quickly (lis per field). Two normally-weighted images 
with no deconvolution were generated for each sub-field (one 
with and one without primary beam correction). The only com- 
putationally demanding steps were generating the images. 

Figure illustrates the results of this process. Only a sim- 
ple model of A has been subtracted so this source is clearly 
visible. Other flux is very likely to be distant sidelobes of A. 
Away from this spurious flux the RMS noise is approximately 
500 jxJy bearrT 1 , approximately double the theoretical prediction 
and similar to that reached by Cornwell et al.. 



5. Discussion 

5.1. The nature of the 1320+299 Complex 

Our observations aimed to clarify the nature of the 1320+299 
complex since the detection (or non detection) of each compo- 
nent on milli-arcsecond scales would further aid their classifica- 
tion. 

Component A was the only one detected by the present ob- 
servations. The image we obtained of component A (Fig. |4]i 
shows the core-jet structure usually found in flat spectrum radio 
quasars. This finding confirms the interpretation suggested by 
Cornwell et al. considering their observations made at lower res- 
olution, in which the A component also shows a core-jet struc- 
ture with the core presenting a flat spectral index. 

The lack of detection of any compact feature in component 
B may confirm the interpretation given by Cornwell et al., who 
have suggested that it is an head-tail galaxy. On the other hand, 
the non-detection of component C, suggested by Cornwell et al. 
to have an appearance similar to that of an edge-brightened hot- 
spot, is still puzzling. 



4 To ensure that the fringe-fit solutions corrected for unmodelled de- 
lays and not any error in the a priori position of component A, the 
astrometric position of component A, phase-referenced to J 1329+3 154 
was carefully measured and another dataset was generated shifted to 
this precise position. 
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(c) (d) 

Fig. 5: Images representing the entire primary beam. Centred on RA 13 h 23 m 2 s .75 Decl. +29°4r32'.'5. Each pixel represents a field 
6.5"x 5". The flux density in Jy beam 1 is shown for the brightest pixel in each image and the RMS of all of the pixels both with and 
without primary beam correction applied: (a) Peak flux density no PB correction; (b) RMS flux density no PB correction; (c) Peak 
flux density PB correction; (d) RMS flux density PB correction. The brightest point (towards the top-right) is the A component, 
attenuated since a simple model of A has been subtracted from the UV data. All flux apart from the brightest pixel in each image is 
spurious flux from the sidelobes of A or noise. 

5.2. Accurate UV shifting 



We have outlined a method whereby a technically inclined as- 
tronomer can image any part of the primary beam of interest. 
Both the accuracy of this method (and the ability of the DiFX 
software correlator to generate visibility datasets with sufficient 
resolution) have been demonstrated for shifts significantly larger 
than those that could be required for Earth-based VLBI with 
parabolic dishes. With the exception of unavoidable instrumen- 
tal effects (the primary beam) the data are not degraded in any 
significant way compared with carrying out separate observa- 
tions. This drastically increases the number of objects which can 
be studied and allows VLBI surveys to be carried out. No corre- 
lation facilities or expertise are required but substantial storage 
and computing resources would be required for larger datasets. 



5.3. Future Work 

The main problem with this approach as it stands is that the 
astronomer still requires access to a sophisticated delay model. 
Code for generating these models, while necessary for correla- 
tion, is not particularly portable nor easy to assimilate into exist- 
ing data analysis software. 

This algorithm has been im plemented with fu ll accuracy 
within the DiFX correlator itself (Deller et al. 20lTb- F° r a large 
subset of wide-field VLBI observations this simplifies the pro- 
cess greatly, obviating the need for large datasets to be written 
to disk at all, and providing the astronomer with a large num- 
ber of standard VLBI datasets. This is particularly appropriate 
for cases where the field of interest has already been studied at 
low resolution with at least equivalent sensitivity to the VLBI 
observation. VLBI surveys of objects detected with other instru- 
ments can provide very useful complementary informations for 
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Fig. 6: An illustration of the non-linearity of the delay (in 
microseconds) with offset for three baselines. Matera-Medicina 
(crosses), Matera-Wettzell (squares), Medicina-Wettzell (cir- 
cles). The top panel shows the delay at offsets between 0.01 and 
100 arcminutes compared to the delay at zero offset. The middle 
panel shows the same delays after the subtraction of the linear 
interpolation between zero offset and the 0.3 arcminute value. 
The bottom panel shows the quotient of the first two plots i.e. 
the fractional error of the linear fit. 



mult-wav elength studies, for exam ple unambiguously identify- 
ing AGN (jMiddeiber^itin|20i3). 

In other cases there will be a benefit to imaging the entire 
primary beam. With the increasing sensitivity of VLBI, it may 
often be the case that pre-existing lower-resolution images are 
not as sensitive as the VLBI studies, at least until instruments 
such as the EVLA and e-MERLIN have completed their legacy 
surveys. Transient phenomena (such as supernovae) may also be 
detected in this way, either with dedicated observations or by 
commensally reprocessing other observations using wide-field 
techniques. Finally, when extended objects are imaged at low 
resolution, large swathes of sky may have to be imaged to find 
compact emission (for example searching for hotspots in jets or 
supernovae and supernova remnants close to a galactic centre. 

For cases where the shifting must be done after correlation, it 
is clear that existing shifting algorithms are inadequate. Figure[6] 
shows the limitation of using a linear delay model (note that the 
linear fit calculated here is identical to that used to calculate the u 
and v coordinates in DiFX). A linear relationship between delay 
and offset (implicitly assumed in interferometry imaging) results 
in an positional error approaching 1 % of the shift within the pri- 
mary beam for a typical VLBI observation. 

As explained in Sect. 12.11 current VLBI datasets contain the 
delay (w) and the partial derivatives of delay with respect to / and 
m (u and v). By adding further derivatives or polynomial func- 
tions of the delay with respect to I and m, the delay across the 
wide field could be fully characterised and accurate UV shifting 
to any point in the wide field could be done at any stage during 
correlation, calibration or imaging. 

At the calibration stage, this three-dimensional (/, m and 
time) model could be refined using in-beam calibrators and/or 
observing techniques which use observations of multiple cali- 
brators. This technique could also be useful for low-frequency 
interferometry where the unmodelled delay is known to vary sig- 
nificantly and non-linearly across the primary beam. 



Calibration using m ultiple weak sourc es across the primary 
beam is described in Garre tt et al ] (12004 We have had some 
early success applying a variation of this technique to other wide 
field datasets, and it is likely to become ever more useful as 
VLBI bandwidths increase. This opens up the possiblity of ob- 
taining noise-limited images and high-accuracy astrometry from 
fields where there is no nearby calibrator. 

In conclusion, advances in correlation and data reduction 
techniques has made VLBI surveys possible and even routine, 
without any loss in sensitivity (or astrometric accuracy) apart 
from that imposed by the instrument. Imaging any or all points 
within the interferometer elements' primary beam with VLBI 
resolution will undoubtedly provide a useful complement to 
wide field surveys at other wavelengths. In addition, data reduc- 
tion techniques developed for wide field VLBI data volumes will 
no doubt be applicable to new wide field and wide band radio 
telescopes. 
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